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Abstract 

The electromagnetic form factor of the pion in the space-like region, and at finite 
temperature, F 7T (Q 2 ,T), is obtained from a QCD Finite Energy Sum Rule. The 
form factor decreases with increasing T, and vanishes at some critical temperature, 
where the pion radius diverges. This divergence may be interpreted as a signal for 
quark deconfinement. 
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The possibility that QCD exhibits chiral-symmetry restoration and quark-gluon decon- 
finement at finite temperature has triggered interest in the thermal behaviour of QCD 
in general, and of hadronic propagators in particular. In this respect, some time ago 
a proposal was made to consider the imaginary part of any hadronic propagator as a 
phenomenological order parameter for the deconfinement phase transition 0. Accord- 
ing to this proposal, with increasing T one should expect resonances to become broader, 
and stable-hadron propagators to develop a non-zero imaginary part. This behaviour 
has been confirmed later for pions, nucleons, and rho-mesons using a variety of theoreti- 
cal approaches 0. Independent phenomenological evidence for the deconfinement phase 
transition in QCD may be obtained e.g. by studying the thermal behaviour of the elec- 
tromagnetic form factor of the pion, F n . In this case one would expect the size of the 
pion to increase with increasing temperature. At the critical temperature the pion radius 
should presumably diverge, indicating quark-gluon deconfinement. A recent calculation 
in the framework of the Nambu-Jona Lasinio model supports this scenario. 

In this note we determine the temperature dependence of F n in the space-like region us- 
ing a Finite Energy QCD Sum Rule (FESR). The pion form factor at T = has been 
extensively studied in the past with FESR, as well as with Laplace transform QCD sum 
rules ||. In order to establish some notation, as well as the T = normalization, we 
briefly describe the method at T = before introducing thermal corrections. 

The appropriate object to study is the three-point function 

IWCP,/,?) = * 2 / d 4 x d A y e^-w) < 0|T(4(x) V x (y) A„(0))|0 > , (1) 

where A^x) = u(x)j^j 5 d(x) is the axial- vector current, V\ is the electromagnetic current, 
and q = p' — p the momentum transfer. On general analyticity grounds, the three-point 
function (1) satisfies the double dispersion relation 

iw^v 2 , q 2 ) = - 2 r ds r ds> ( p z x{ *;?;®\ , (2) 

TV 2 JO JO [S + p Z )(s' + p' 2 ) 

defined up to subtractions, which are disposed of by Laplace improving the Hilbert trans- 
form, or by considering FESR. The QCD sum rule program at T = || may be phrased 
as follows. The left hand side of (2) in the Euclidean region can be calculated in QCD 
through the Operator Product Expansion (OPE), to any desired order in perturbation 
theory, and parametrizing non-perturbative effects in terms of vacuum expectation val- 
ues of the same quark and gluon fields entering the QCD Lagrangian. Contact with the 
hadronic world is made by writing the spectral function appearing in the right hand side 
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of (2) in terms of all possible hadronic intermediate states contributing to the three-point 
function. Normally, a reasonable saturation of the dispersion relation is achieved by the 
ground state, followed after some threshold energy s by a hadronic continuum modelled 
by perturbative QCD. 



The correlator (1) involves quite a few structure functions, associated with all the Lorentz 
structures that can be formed with the available four-momenta. In principle, it should not 
matter which particular structure one chooses to project the pion form factor. We follow 
|| in choosing the combination P^P u Px, where P = p + p'. In this case, the hadronic 
spectral function in the chiral-limit reads 

p(s, s', Q 2 )\had = \ f 2 F 7T (Q 2 )5(s)5(s') + p(s, s', Q 2 )\ QC d[1 ~ 0(s - s - s')} , (3) 



where s signals the onset of the continuum, /„. ~ 93 MeV, and 

P(S, S', Q 2 )\ QCD = JL [ 3A(a . + Q2 )(x + 2Q 2 ) _ A 2 _ 5Q 2 {x + g2) 3] ? (4) 

to one- loop order (and in the chiral-limit), with 

A = y 2 + Q 2 (2x + Q 2 ) , (5) 



and x = s + s', y = s — s'. Since we are interested in writing the lowest moment FESR 
for F n , i.e. 



] i-sq rx 

Fn(Q 2 ) = / dx dy p(x,y,Q 2 )\ QC D, (6) 

rather than a Laplace transform QCD sum rule, the non-perturbative power corrections 
entering the OPE are of no concern to us here (they contribute to higher moment FESR). 
The integration region in (6) has been chosen to be a triangle in the (s,s') plane, with 
base and height equal to sq. Other choices of the integration region, e.g. a square region 
of side s\ ~ so/v^2, give similar results. The solution to the FESR (6) is 

F ^ Q2) = 167r 2 /2(l + Q 2 °/2 So ) 2 • (7) 



Although not evident from (7), it is important to realize that this analysis is only valid in 
the region Q 2 > 1 GeV 2 , where one expects a reasonable convergence of the OPE. This 
limitation is of no relevance if one is only interested in the thermal behaviour of the ratio 
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F n (Q 2 ,T)/F n (Q 2 ,0). In any case, as shown in J|, Eq.(7) provides a reasonable fit to the 
experimental data in the region Q 2 ~ 1 — 4 GeV 2 , if so — 1 GeV 2 . 

We now proceed to derive the thermal corrections to (6). The vacuum average in Eq.(l) 
is to be replaced by the Gibbs average 

« A ■ B ■ • • »= exp(-£„/T) < n\ A ■ B • ■ ■ \n > /Tr(exp(-H/T)) , (8) 

n 

where |n > is any complete set of eigenstates of the (QCD) Hamiltonian, e.g. hadronic 
states, quark-gluon basis, etc.. If one is to extend smoothly the T = QCD sum rule 
program to finite temperature, then the natural choice for \n > should be the quark-gluon 
basis, as first proposed in 0. In fact, assuming the validity of the OPE at T ^ and 
invoking QCD-hadron duality at T = 0, raising the temperature by an arbitrary small 
amount would induce very little change in the hadronic spectrum, and in the expecta- 
tion values of the QCD operators in the OPE. Hence, it is reasonable to expect that the 
inter-relationship between QCD and hadronic parameters effected by duality will remain 
valid. An abrupt dissapearance of this inter-relationship as soon asT ^ appears highly 
unlikely. As to the validity of the OPE at finite temperature, one is faced with the same 
problem as when T = 0, i.e. no rigorous proof can be given since one is not able to solve 
QCD analytically and exactly. Instead, other field theory models which are exactly solv- 
able have been used to argue for the validity of the OPE at T = J7|. This analysis has 
been extended recently to T / ||, and it shows the same level of supportive evidence 
as in the T = case. 



The QCD program at T ^ as outlined above leads to some interesting results. For 
instance, considering the two-point function associated with the axial-vector current, a 
lowest moment FESR relates sq(T) to f n (T) in such a way that sq(T) vanishes at some 
temperature T d @-[10|. This may be interpreted as the critical temperature for deconfine- 



ment. In fact, with increasing T, and as resonances begin to melt, the hadronic spectrum 
should smooth out and s (T) should decrease. Close to T d the spectral function is then 
described almost entirely by the quark-gluon degrees of freedom at all energies. Depend- 
ing on the expression for / W (T) used as input, T d could be smaller or almost the same 
as the critical temperature for chiral-symmetry restoration, T c . For instance, using chiral 
perturbation theory for f n (T) leads to T d < T c ||, but this input is only valid at low 
temperatures. A more refined analysis gives |10| 



s (T) U{T) 
\ *o(0) W) 



(9) 



and Trf only slightly lower than T c (given the accuracy of the method this difference is of 
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not much significance). 



We have calculated the spectral function (4) at finite temperature using the Dolan-Jackiw 
formalism. After substitution in (16), the result can be expressed as 

1 pso( T ) r x 

F *(Q 2 iT) = —— / dx / dyp(x,y,Q 2 )\ QC DF(x,y,Q 2 ,T), (10) 
f4(T) Jo J-x 

with 

F(x, y, Q 2 , T) = 1 - ni - n 2 - n 3 + n x n 2 + + n 2 n 3 , (11) 



"j "2 = I I^W^~y^' I ■ (12) 



, , Q 2 + (x-y)/2 , , 

n 3 = ^F 7= , (13) 

and rip is the Fermi thermal factor. In the equations above we have chosen a frame such 
that p M = (no, 0), and p'^ = (a/, p'), in which case 



w = J^ , o/ = ^±£. (14) 
V 2 2 fZ& 



We have explicitly checked that the ratio 



R{T) - TO (15) 



is essentially insensitive to other choices of frames. For instance, one may choose = 
((j, p), and p'^ = (u/, — p), which leads to different arguments in the thermal factors, but 
roughly the same ratio R(T). The temperature behaviour of the ratio (15) is shown in 
Fig.l for Q 2 = 1 GeV 2 (solid curve), and Q 2 = 3 GeV 2 (broken curve), where we have used 



(9) together with the results of JlO for f n (T). This result for R(T) is in nice agreement 
with the expectation that as the temperature increases, F n should decrease and eventually 
vanish at the critical temperature for deconfinement T<j. 
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Although the OPE breaks down at small values of Q 2 , one may still extrapolate the ratio 
(15) into this region just to study the qualitative temperature behaviour of the root- 
mean-square radius ratio < r 2 . > T / < r\ > . Doing this, we find that this ratio increases 
monotonically with T, doubling at T/Tj ~ 0.8, and diverging at the critical temperature. 
This divergence of < r\ may be interpreted as a signal for quark deconfinement. In 
fact, the behaviour of < r\ >t can be traced back to the temperature behaviour of the 
asymptotic freedom threshold sq(T). As sq(T) decreases with increasing T, a signature 
of quark deconfinement, the root-mean-square radius of the pion increases. This result is 
in qualitative agreement with that obtained in the framework of the Nambu-Jona Lasinio 
model § 
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Figure Captions 

Figure 1: Temperature beahviour of the ratio (15) for Q 2 = 1 GeV 2 (solid curve), and 
Q 2 = 3 GeV 2 (broken curve). 
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This figure "figl-l.png" is available in "png" format from: 
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